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ABSTRACT 

The T Tauri star PTFO 8-8695 exhibits periodic fading events that have been interpreted as the transits of 
a giant planet on a precessing orbit. Here we present three tests of the planet hypothesis. First, we sought 
evidence for the secular changes in light-curve morphology that are predicted to be a consequence of orbital 
precession. We observed 28 fading events spread over several years, and did not see the expected changes. 
Instead we found that the fading events are not strictly periodic. Second, we attempted to detect the planet’s 
radiation, based on infrared observations spanning the predicted times of occultations. We ruled out a signal 
of the expected amplitude. Third, we attempted to detect the Rossiter-McLaughlin effect by performing high- 
resolution spectroscopy throughout a fading event. No effect was seen at the expected level, ruling out most 
(but not all) possible orientations for the hypothetical planetary orbit. Our spectroscopy also revealed strong, 
time-variable, high-velocity H« and Ca H & K emission features. All these observations cast doubt on the 
planetary hypothesis, and suggest instead that the fading events represent starspots, eclipses by circumstellar 
dust, or occultations of an accretion hotspot. 

Subject headings: planetary systems - stars: individual (PTFO 8-8695) - stars: pre-main sequence 


1. INTRODUCTION 

The discovery of close-in giant planets around very young 
stars — less than a few million years old — would provide 
precious information about the timing of planet formation, the 
structure of newborn planets still cooling and contracting, and 
the mechanism for shrinking planetary orbits and creating hot 
Jupiters. Currently the only candidate for such an object is 
PTFO 8-8695b, found by van Eyken et al. pOT2]) (hereafter 
VE+12). PTFO 8-8695 is a T Tauri star in the Orion-OBla 
region, with a mass of ^0.4 M 0 , a spectral type of M3, and 
an estimated age of 3 Myr (Briceno et al. |2005| ). In addition to 
the quasi-sinusoidal variability characteristic of T Tauri stars, 
this star was found to exhibit periodic fading events, during 
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which the star dims by a few percent for an interval of about 
1.8 hours. VE+12 reported on these and other properties of 
the system, and advanced the hypothesis that the fading events 
are transits of a close-in giant planet. 

However, the planetary interpretation is not secure. The 
system has some puzzling properties that seem incompatible 
with the planet hypothesis, or at least demand that the sys¬ 
tem has somewhat exotic properties. In the first place, the 
“transit” light curves do not have the customary morphology. 
They were seen to vary in depth and duration over a timespan 
of a year, and in some cases to lack the expected symm etry 
around the time of minimum light. Barnes et al. ( f2Q13| ) pro¬ 
posed that these changes are caused by a large misalignment 
between the planet’s orbit and the star’s equatorial plane. This 
misalignment, when combined with an asymmetric intensity 
profile on the stellar disk due to gravity darkening, can pro¬ 
duce asymmetric transit light curves. Furthermore, the mis¬ 
alignment leads to nodal precession of the orbit, which could 
explain the secular changes in morphology. 

By itself this would not be too unusual. Hot Jupiters with 
spin-orbit misalignments are now commonplace (see, e.g. Al¬ 
brecht et al. |2012|), and nodal precession has been observed 


in at least one other misaligned system (Szabo et al. [201 1 
|2012| . Barnes et al. ( |2013| constructed a model that quan 
titatively fits the two light curves measured by VE+12 in 
2009 and 2010. However, in the case of PTFO 8-8695 the 
“transit” period is equal to the stellar rotation period (as esti¬ 
mated from the quasi-sinusoidal variability): both are consis¬ 
tent with 0.448 days or 10.8 hours. It seems strange that the 
system would have reached spin-orbit synchronization with- 


vnc J 

out also achieving spin-orbit alignment]^] The coincidence 


15 Kamiaka et al. <2015] explored models in which the orbital and rotation 
periods are not necessarily synchronized, under the premise that the stellar 
rotation period could have any value up to 16 hours (an upper limit set by the 
measured vsinL). However, the quasi-sinusoidal flux variations outside the 
fading events are likely due to rotation, and have a period that agrees with 
that of the “transit” events to within a percent. Thus it seems unnecessary to 
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between the “transit” and rotation periods raises the possi¬ 
bility that the fading events are actually due to starspots, or 
eclipses by a corotating structure within a circumstellar disk 
or accretion flow. 

Another striking property of PTFO 8-8695 is that the plan¬ 
etary radius inferred by VE+12 was 1.9 Rj m making it essen¬ 
tially tied with WASP-17b (Triaud et al. |2010| for the largest 
known planetary radius. Perhaps this should be expected, for 
a planet that is still contracting from an initially distended 
state. Somewhat more worrying is that the orbital period of 
10.8 hours is within or at least near the Roche limit for a gas 
giant (Rappaport et al. |2013| ). This suggests that the planet 
would be actively losing mass through Roche lobe overflow. 

The conventional way to confirm the existence of a transit¬ 
ing planet is to detect the expected radial-velocity variation 
of the host star. VE+12 attempted to detect such a signal but 
were foiled by the spurious radial-velocity variations caused 
by stellar activity, which are larger than the amplitude of the 
expected orbital velocity. Even if an apparently sinusoidal 
radial-velocity signal were detected, it would be difficult to 
ascertain whether the signal is planetary in origin or arises 
from stellar activity, due to the coincidence between the “tran¬ 
sit” and rotation periods. 

Given the high scientific stakes, we attempted three less 
conventional tests of the planetary hypothesis: 

1. According to the gravity-darkening model of Barnes 
et al. ( [2013] ), continued nodal precession should pro¬ 
duce variations in the asymmetry, duration, and depth 
of the fading events, with a period of a few years. 
Barnes et al. ( |2013| ) also predicted that there should 
be intervals of several months during which the fading 
events cease, because the planet’s trajectory does not 
cross the face of the star. Therefore, we undertook time- 
series photometry of as many fading events as possible 
over a timespan of several years, to detect the expected 
changes in morphology. 

2. Close-in giant planets emit relatively strongly at in¬ 
frared wavelengths, due to a combination of reflected 
starlight and the planet’s own thermal radiation. There¬ 
fore, we attempted to detect the loss of light when the 
planet is hidden by the star, by performing time-series 
infrared photometry spanning the expected times of oc- 
cultations (halfway between transits). 

3. A key premise of the planet hypothesis is that the or¬ 
bit is misaligned with the stellar equator. In con¬ 
trast, starspots move in a prograde direction, aligned 
with stellar rotation. The angle between the trajec¬ 
tory of a transiting feature and the (sky-projected) stel¬ 
lar equator can be measured by observing the Rossiter- 
McLaughlin effect (Rossiter |1924[ McLaughlin 11924| . 
Therefore, we undertook high-resolution optical spec¬ 
troscopy throughout a fading event to measure the 
spin-orbit angle of whatever is apparently blocking the 
starlight. We also used the spectra to check for time 
variations in the sky-projected rotation rate (vsinz*), 
which would be expected if the star is precessing. 

This paper is organized as follows. Section[2]presents time- 
series photometry of the candidate transits and occultations, 
using several ground-based telescopes and an archival obser¬ 
vation with the Spitzer Space Telescope. Section [3] presents 

consider non-synchronized models. 



Year of observation 

Fig. 1.— Dates of observations of fading events with the FLWO 1.2m and 
TRAPPIST 0.6m telescopes, along with previous observations by VE+12. 
Also indicated are the dates of the observations with Spitzer, Magellan, and 
Keck. 


our time-series spectroscopy and our attempt to detect the 
Rossiter-McLaughlin effect. Section |4] analyzes the preced¬ 
ing results and their implications for the planetary hypothesis 
as well as other possible explanations for the fading events. 


2. TIME-SERIES PHOTOMETRY 


2.1. Overview 

We conducted time-series photometric observation s of fad- 
ing events between 2012 and 2015. Below, in § |2.2| we 
present ground-based observations of 26 different events. Fig¬ 
ure [I] shows their distribution in time, and Figure [2] shows the 
light curves (including the 7 highest-quality light curves pre¬ 
sented previously by VE+12, for reference). In a few cases we 
observed the event through multiple broadband filters. The re¬ 
sulting multi-band light curves are shown in Figure[6] We also 
observed a candidate occultation at infrared wavelengths with 
one of the Magellan 6.5m telescopes; those data are described 
in § 2.3. 1| and plotted in Figure |7] Finally, we analyzed the 
available Spitzer data, spanning a fading event as well as an 
expected occultation. Those data are described in § |2.3.2| and 
shown in Figure [9] The dates of the Spitzer and Magellan ob¬ 
servations are also indicated on Figure [T] along with the Keck 
spectroscopic observations described in Section [3] 


2.2. Ground-based observations of fading events 

We observed 13 fading events with the 1.2m telescope at 
the Fred Lawrence Whipple Observatory (FLWO) on Mt. 
Hopkins, Arizona. The instrument, Keplercam, has a single 
4096 x 4096 CCD with a 23"1 field of view. All the events 
were observed through an i' filter. For the events between 
2014 Jan 9-18, we interleaved the i' - band exposures with g f - 
band exposures, although the g'-b and data were only useful in 
two cases. Calibration was performed using standard IRAff^] 
procedures, including bias and flat-field corrections. The time 
stamps were placed on the BJD T db system using the code by 


16 The Image Reduction and Analysis Facility (IRAF) is distributed by the 
National Optical Astronomy Observatory, which is operated by the Associa¬ 
tion of Universities for Research in Astronomy (AURA) under a cooperative 
agreement with the National Science Foundation. 
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TABLE 1 

Best-FITTING PARAMETERS of fading events (see Eq.[TJ. 


UT date 

Filter 

^min 

No. of 
data points 

Time of minimum light, to 
[HJD-2,455,200] 

Fractional loss 
of light, 5 

Duration, w 
[days] 

2010 Jan 5 

R 

193.4 

228 

1.8053 ±0.0014 

0.0442 ±0.0052 

0.0282 ±0.0033 

2010 Jan 6 

R 

216.5 

244 

2.6961 ±0.0017 

0.0789 ±0.0208 

0.0349 ±0.0059 

2010 Jan 9 

R 

217.0 

191 

5.8431 ±0.0014 

0.0979 ±0.0142 

0.0370 ±0.0042 

2010 Dec 9 

R 

169.9 

270 

339.9069 ±0.0006 

0.0498 ±0.0023 

0.0190 ±0.0013 

2010 Dec 10 

R 

248.4 

298 

340.8020 ±0.0006 

0.0425 ±0.0011 

0.0167 ±0.0007 

2010 Dec 13 

R 

191.1 

271 

343.9467 ±0.0010 

0.0558 ±0.0038 

0.0229 ±0.0017 

2010 Dec 14 

R 

473.2 

288 

344.8401 ±0.0009 

0.0412 ±0.0018 

0.0158 ±0.0012 

2012 Sep 4 

I+z 

194.0 

209 

974.8427 ±0.0005 

0.0510 ±0.0027 

0.0160±0.0011 

2012 Sep 13 

I+z 

388.8 

265 

983.8121 ±0.0004 

0.0521 ±0.0014 

0.0138 ±0.0006 

2012 Nov 23 

I+z 

240.1 

243 

1054.6614 ±0.0009 

0.0189 ±0.0013 

0.0128 ±0.0014 

2012 Dec 11 

I+z 

94.2 

94 

1072.5994 ±0.0009 

0.0366 ±0.0024 

0.0132 ±0.0015 

2012 Dec 14 

I+z 

154.0 

176 

1075.7379 ±0.0006 

0.0298 ±0.0014 

0.0098 ±0.0008 

2012 Dec 15 

I+z 

152.8 

182 

1076.6357 ±0.0005 

0.0338 ±0.0029 

0.0128 ±0.0013 

2013 Oct 12 

i' 

9.7 

54 

1377.9469 ±0.0009 

0.0324 ±0.0024 

0.0099 ±0.0012 

2013 Nov 23 

I+z 

102.7 

93 

1419.6524 ±0.0022 

0.0177 ±0.0029 

0.0114 ±0.0030 

2013 Dec 1 

I+z 

245.9 

195 

1427.7243 ±0.0007 

0.0316 ±0.0026 

0.0163 ±0.0016 

2013 Dec 11 

I+z 

190.3 

233 

1437.5864 ±0.0007 

0.0385 ±0.0030 

0.0153 ±0.0016 

2013 Dec 14 

I+z 

255.1 

262 

1440.7266 ±0.0005 

0.0344 ±0.0016 

0.0139±0.0011 

2014 Jan 9 

i f 

127.6 

46 

1466.7304 ±0.0007 

0.0385 ±0.0016 

0.0156 ±0.0012 

2014 Jan 14 

i' 

183.0 

57 

1471.6634 ±0.0011 

0.0337 ±0.0025 

0.0128 ±0.0015 

2014 Jan 17 

i' 

245.0 

59 

1474.8039 ±0.0010 

0.0377 ±0.0042 

0.0160 ±0.0023 

2014 Jan 18 

i' 

177.0 

57 

1475.7000 ±0.0007 

0.0375 ±0.0014 

0.0144 ±0.0009 

2014 Jan 19 

I+z 

253.6 

240 

1476.5974 ±0.0007 

0.0326 ±0.0023 

0.0167 ±0.0016 

2014 Jan 23 

i' 

331.9 

136 

1480.6289 ±0.0004 

0.0628 ±0.0044 

0.0218 ±0.0013 

2014 Jan 23 

I+z 

549.4 

356 

1480.6339 ±0.0010 

0.0295 ±0.0018 

0.0149 ±0.0015 

2014 Feb 5 

i' 

404.4 

111 

1493.6353 ±0.0005 

0.0349 ±0.0013 

0.0140 ±0.0008 

2014 Feb 9 

i' 

330.9 

107 

1497.6729 ±0.0005 

0.0402 ±0.0013 

0.0132 ±0.0008 

2014 Feb 13 

i' 

246.4 

99 

1501.7075 ±0.0005 

0.0451 ±0.0022 

0.0147 ±0.0010 

2014 Feb 22 

i' 

476.4 

83 

1510.6778 ±0.0007 

0.0463 ±0.0027 

0.0148 ±0.0013 

2014 Nov 17 

i' 

208.1 

119 

1778.7654 ±0.0007 

0.0304 ±0.0020 

0.0116±0.0012 

2014 Nov 29 

i' 

259.0 

145 

1790.8761 ±0.0009 

0.0185 ±0.0014 

0.0088 ±0.0011 

2014 Dec 27 

i' 

199.7 

118 

1818.6748 ±0.0007 

0.0222 ±0.0012 

0.0094 ±0.0008 

2015 Feb 27 

I+z 

200.0 

206 

1880.5567 ±0.0009 

0.0206 ±0.0017 

0.0087 ±0.0012 


TABLE 2 

Best-fitting parameters of phase-folded light curves (see Eq.[2) 


Light curve 

Filter 

Time of minimum light, to 
[days] 

Fractional loss 
of light, S 

Ingress duration, w\ 
[days] 

Egress duration, W 2 
[days] 

FLWO average 
TRAPPIST average 

i' 

I+z 

0.0029 ±0.0028 
0.0018 ±0.0018 

0.0360 ±0.0017 
0.0336 ±0.0010 

0.0166 ±0.0023 
0.0149 ±0.0014 

0.0124 ±0.0018 
0.0128 ±0.0013 


Eastman et al. ( |2010| ). Circular aperture photometry was per¬ 
formed with the Interactive Data Language (IDL). 

Another 13 events were observed with the 0.6m TRAnsit- 
ing Planets and Planeteslmals Small Telescope (TRAPPIST), 
located at ESO’s La Silla Observatory in Chile. This telescope 
is equipped with a thermoelectrically-cooled 2048 x 2048 
CCD with a 22' field of view (Gillon et al. 2011 ; Jehin et al. 
[2011| ). The observations were conducted with a custom “7+z” 
filter, which has transmittance >90% between 750-1100 nm. 
We refer the reader to Gillon et al. ( |2013| ) for descriptions 
of the procedures for observing and data reduction. The two 
events of 2012 Dec 14 and 15 were simultaneously observed 
with a Gunn r' filter, using EulerCam on the 1.2m Euler- 
Swiss Telescope at the La Silla site, Chile. EulerCam uses 
a 2048 x 2048 CCD with a field of view of 14"7 on a side. 
For details on the instrument and data reduction procedures, 
please refer to Lendl et al. ( 2012| ). 


A single event on 2014 Jan 19 was observed with the 6.5m 
Magellan I (Baade) telescope at Las Campanas Observatory 
in Chile. The same event was observed simultaneously by 
TRAPPIST in the I + z band. With Magellan, we observed in 
the H band using FourStar, a 2048 x 2048 infrared array with 
a 1(T9 square field of view. The data were reduced with IRAF 
and IDL procedures similar to those used on the FLWO data. 

In all cases the flux of PTFO 8-8695 was divided by the 
summed flux from several reference stars, leading to the light 
curves plotted in Figure [2] This figure also shows the 7 light 
curves presented by VE+12 that cover the entire fading event; 
those observations were performed with the 1.2m Palomar 
telescope and an R filter. 

Outside of the fading events, the star varies gradually by 
~0.1 mag over several hours, in a manner consistent with its 
young age and late spectral type. Superimposed on those rela¬ 
tively gradual variations are periodic transit-like fading events 
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Phase 

Fig. 2.— Time-series photometry of fading events of PTFO 8-8695. The vertical scale is the same for all events; vertical offsets have been applied to separate 
the different time series. Included in this plot are the 7 highest-quality light curves from VE+12 (red), as well as 26 new light curves from the FLWO 1.2m 
telescope (blue) and the TRAPPIST 0.6m telescope (magenta). The FLWO observations between 2014 Jan 9-18 were conducted in both i' and g' band. 
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lasting no more than about 2 hours. The depth and duration 
of the fading seems to vary from event to event. To derive the 
basic phenomenological parameters of the dimming events — 
depth, duration, and time of minimum light — we fitted a pa¬ 
rameterized model describing both the gradual out-of-transit 
variations as well as the transit-like loss of light. We modeled 
the gradual variations as a polynomial function of time (2nd 
or 3rd order, depending on the event). The additional loss 
of light during the fading event, relative to the polynomial- 
corrected out-of-transit flux, was modeled as 


A fit) = 5 sech 


. w . 


25 

e (t-to)/w _|_ e -(t-to)/w 5 


( 1 ) 


where 5 is the maximum fractional loss of light (the “transit 
depth”), w is the duration, and is the time of minimum light. 
We chose this model instead of a more physically-motivated 
transit model (e.g., the model presented by Mandel & Agol 
|2002| ) because the asymmetries and other odd features in the 
light curves do not fit the standard models. Hence there is no 
advantage in fitting the physical model when a much simpler 
model can provide estimates of the basic transit parameters. 
One might be able to fit the data with a model based on tran¬ 
sits of an oblate, oblique, precessing, gravity-darkened star 
(Barnes [2009 ), but such a model is far more demanding com¬ 
putationally. Our analytic model suffices to estimate the basic 
parameters of each event. 

Figured] shows the light curves after dividing out the best¬ 
fitting polynomial functions. This gives a clearer view of the 
“transits” with most of the long-term trends removed. Ta¬ 
ble jllgives the model parameters, as well as the value of x^in 
ancTthe number of data points in each time series. In most 
cases, Xmin i s t0 ° large to be statistically acceptable, i.e., the 
simplified model of Eq. [I] does not fit the data to within the 
photometric uncertainties. For this reason, Table [T] does not 
report the formal parameter uncertainties defined by the usual 
criterion A% 2 = 1. Rather, the reported parameter uncertain¬ 
ties have been enlarged by a factor of \J X^in/^of where N d 0 f 

is the number of degrees of freedom. These enlarged uncer¬ 
tainties were also adopted for our subsequent calculations. 

Figure [4] shows the measured depths and durations. When 
the transits are deeper, they also tend to have longer durations; 
the measured depth and duration have a Pearson correlation 
coefficient of 0.84 (p < 10 -5 ). 

The VE+12 light curves showed strong asymmetries in at 
least 5 out of the 7 complete light curves. None of our new 
light curves show strong asymmetries, at least not as clearly 
as was seen by VE+12. We tried fitting a model in which 
the rates of brightness variation are not symmetric about the 
time of minimum light, by using a non-standard variant of the 
hyperbolic secant function (see, e.g., Ruan et al. |2000| ): 


a/=-^ 


25 


e (t-to)/w] _|_ e -(t-to)/w 2 ' 


( 2 ) 


The asymmetric model does not seem to improve the quality 
of the fit to a significant degree. The number of cases for 
which the fitted asymmetry obeyed w\ > 1+2 (more prolonged 
“ingress”) was nearly the same as the number of cases with 
W 2 > i+i, without any obvious pattern. Figure [5] shows the 
FLWO and TRAPPIST data as a function of t-tp (coverted to 
hours) after some averaging in time to increase the signal-to- 
noise ratio. In these averaged light curves there does seem to 
be a slight asymmetry, with a longer “ingress” than “egress” 


in both cases. The best-fitting asymmetric model is shown in 
Figure [5] and the parameters are given in Table [2] 

Another finding is that the loss of light is usually strongly 
chromatic, as shown in Figurejb] For those two cases in which 
we observed the same event in both the i' and g' bands, we 
found the loss of light to be 30-50% larger in the g' band. For 
the single case in which we observed in both 77 and 7 + z, the 
loss of light was also ^-40% larger in the bluer band. We also 
observed two events simultaneously in the Gunn r' and 7 + z 
bands. On the first night the loss of light was 20-30% larger 
in the bluer band. On the second night, the loss of light in r' 
was essentially the same as in 7+z. 

2.3. Ground-based observations of predicted occultations 

The infrared emission from close-in giant planets is often 
strong enough that it is possible to detect the decrement in to¬ 
tal flux when the planet is hidden by the star. The detection of 
such an occultation would be strong evidence supporting the 
planet hypothesis. We searched for the predicted occultations 
in the Magellan and Spitzer time-series infrared photometry, 
as discussed below. 


2.3.1. Magellan observations 

The 6.5m Magellan I (Baade) telescope and FourStar were 
also used to monitor PTFO 8-8695 for approximately 6 hours 
on 2014 Jan 21, spanning the predicted time of a planetary 
occultation. The prediction was based on the assumption that 
the planet’s orbit is circular, i.e., that the occultations occur 
exactly halfway between transits. The observations were con¬ 
ducted in the 77 ban d an d the data were processed in the man¬ 
ner described in § |2.2| Figure [7] shows the resulting light 
curve. No occultation signal is evident, although some grad¬ 
ual variability is seen. 

We determined an upper limit on the relative brightness 
of the planetary dayside by fitting an occultation model 
to the data. The model had the same total duration and 
ingress/egress durations as the “transit” light curve observed 
with Magellan two days earlier. As with the transits, we fitted 
the gradual variations with a polynomial function of time (in 
this case, a cubic function). We used a Markov Chain Monte 
Carlo (MCMC) algorithm to calculate the posterior probabil¬ 
ity distribution for the parameters describing the cubic func¬ 
tion as well as the loss of light during the occultation. The 
minimum x 2 value was 2082.5, with 1876 degrees of free¬ 
dom, indicating a statistically unacceptable fit; the cubic func¬ 
tion is evidently not a completely satisfactory description of 
the observed flux variations. Rather than develop more elabo¬ 
rate models we simply inflated the parameter uncertainties by 

the factor yjx min/^dof- The resulting occultation depth was 

5qqq = 0.00024 ±0.00016, corresponding to a 3a upper limit 
of 5 occ < 0.00072. 

To decide if this upper limit rules out the planetary hypoth¬ 
esis, we need to know the expected occultation signal. Fol¬ 
lowing the usual simplified model for transiting planets, the 
fractional loss of light during an occultation is 


^OCC — Ag 



2 1 (R P \ 2 .li B ^ dX 

\ R J ± B\(T t )<JX 


The first term is due to reflected starlight, in which A g is the 
geometric albedo, and R p and R* are the planetary and stellar 
radii. The second term is due to the planet’s thermal emis¬ 
sion, in which the observing bandpass extends from Ai to A 2 
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Fig. 3.— Same as Figure [2] but after “flattening” the out-of-transit flux variation by dividingby the best-fitting polynomial function of time. Flux values are 
normalized to unity outside of the transit. The black curves are the best-fitting models (see Eqn.Q], from which we derived the transit times, depths, and durations 
that are reported in Table[l] 
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Fig. 4.— Depths and durations of the fading events, estimated for each in¬ 
dividual event by fitting a simple analytic model (Eq.[2j. There is a positive 
correlation between depth and duration. The best-fit straight line is shown as 
a dashed line. The Pearson correlation coefficient and its statistical signifi¬ 
cance are given in the bottom right corner. 



t-to (hr) 

Fig. 5.— Top. —Phase-folded light curve based on data from all the fading 
events observed with the FLWO 1.2m telescope in the i' band. Bottom .— 
Same, but for the I+z data obtained with the 0.6m TRAPPIST telescope. In 
both cases the data were placed into 50 time bins spanning a 4-hour period 
bracketing the expected transit time. The grey line represents the best-fit 
asymmetric model for each light curve. The error bars represent the standard 
deviation of the mean in each time bin. 

(1.48-1.76 /im, for H band), B\(T) is the Planck function, T\ 
is the star’s effective temperature, and T p is the planet’s day- 


side effective temperature. The latter is calculated from the 
condition of radiative equilibrium, giving 



where A is the Bond albedo and / is a dimensionless number 
depending on the manner of radiation. If the entire surface 
radiates isotropically as a blackbody, then f -4. If instead the 
dayside radiates uniformly and the nightside radiation can be 
neglected, then f - 2. Furthermore, if the angular dependence 
of the planet’s radiation is assumed to follow Lambert’s law, 
then A g = 2A/3. 

In this case it is difficult to establish the key parameters 
(. R p /R*) 2 and R*/a, because of the changing morphology of 
the fading events with time and wavelength. Barnes et al. 
(2013) found (Rp/R*) 2 « 0.027 and R*/a « 0.58 using a 
model incorporating the effects of gravity darkening and or¬ 
bital precession. Using those parameters, we calculate the 
expected value of £ occ and plot it as a function of the Bond 
albedo in Figure [8j for both / = 4 and f - 2. The expected 
£occ ranges from a minimum value of 0.0028 for / = 4 and 
A — 0, to a maximum value of 0.0061 for / = 2 and A — 1. 
Such large occultation depths are ruled out by our Magellan 
observations. 

2.3.2. Spitzer observations 

The Spitzer Space Telescope monitored PTFO 8-8695 on 
2012 Apr 28 for about 12 hours, slightly longer than a full 
photometric period. The data consist of 1369 full-array im¬ 
ages from the IRAC detector operating at 4.5 pm, with an 
integration time of 26.8 seconds. The data were calibrated 
by the Spitzer pipeline version S19.1.0. These observations 
were carried out in non-cryogenic mode under program no. 
80257 (PI: Stauffer), and are publicly available on the Spitzer 
Heritage Archive databasep] 

We converted the data from the Spitzer units of specific 
intensity (MJy sr -1 ) into photon counts, and then performed 
IRAF aperture photometry on each subarray image. Best re¬ 
sults were obtained with an aperture radius of 2.5 pixels and 
a background annulus extending from 11 to 15.5 pixels from 
the center of the point-spread function (PSF). The center of 
the PSF was measured by fitting a two-dimensional Gaussian 
function to each image. At this stage, 30 discrepant fluxes 
were discarded by applying a 5 a median clipping algorithm. 

Next we needed to remove the apparent flux variations as¬ 
sociated with motion of the image on the detector, the main 
source of systematic effects in time-series photometry with 
the IRAC InSb arrays (Knutson et al. |2008| ). This effect is 
caused by the combination of (i) the coarse sampling of the 
PSF, (ii) the significant inhomogeneity of the pixels, and (iii) 
fluctuations in the telescope pointing. To mitigate this ef¬ 
fect we chose the Bi-Linearly-Interpolated Sub-pixel Sensi¬ 
tivity (BLISS) mapping method presented by Stevenson et al. 
( |2012| ). In this method, the flux data themselves are used to 
constrain a model for the subpixel sensitivity variations. In 
our implementation of the method, the detector area probed 
by the PSF center was divided into a 13 x 13 grid. With this 
degree of sampling, the PSF center visited each grid point at 
least 10 times throughout the course of the observations. 

17 While this manuscript was in preparation, Ciardi et al. (20151 reported 
an independent analysis of these same data. 
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Fig. 6.— Multiband observations of fading events. Top row .—Simultaneous observations in the r' and I+z bands (open blue and filled red circles, respectively). 
Middle row .—Interleaved observations in the g' and i' bands (open blue and filled red circles, respectively). Bottom .—Simultaneous observations in th e I+z and 
H bands (open blue and filled red circles, respectively). In all but one case, the loss of light is greater in the bluer bandpass. 


Figure [9] shows the Spitzer light curve after BLISS correc¬ 
tion. The variability can be described as the combination of 
quasi-sinusoidal variation with a period of «0.5 day, and the 
transit-like dip in brightness at the expected time, with an am¬ 
plitude of 0.5% and a duration of approximately 1.4 hours. 
No occultation is seen at the expected time (0.224 days after 
the transit). Figure [TO] gives a better view of the transit-like 
event, and the data surrounding the predicted time of occul¬ 
tation. In this figure, only the data within 0.1 days of each 
event are shown, and the data have been rectified by fitting 
a quadratic function of time to data outside of the event and 
then dividing by the best-fitting function. 

Just as with the Magellan light curve, we determined an 
upper limit on the occultation loss of light by fitting a param¬ 
eterized model to the data. The model included a quadratic 
function of time to describe the out-of-occultation variations. 
The occultation model was required to have the same dura¬ 
tions between first, second, third and fourth contacts as ob¬ 
served earlier with Magellan. The loss of light S occ was a free 
parameter. We used an MCMC algorithm to calculate the pos¬ 
terior probability distribution of S occ and the parameters of the 
quadratic function. The minimum x 2 value was 565.9, with 
535 degrees of freedom. The result for the occultation depth 
was £occ = -0.0008 dz 0.0006, i.e., the best-fitting value cor¬ 


responds to a brightness increase rather than a loss of light. 
This corresponds to a 3cr upper limit of S occ < 0.0010. Again, 
as illustrated in Figure [8] the upper bound on S occ given by 
Spitzer is smaller than the occultation depth implied by the 
parameters of the Barnes et al. (2013 ) model. 

Careful inspection of Figure 9] sEows a candidate flux dip 
of centered around a time coorainate of 0.425, with an am¬ 
plitude of ^0.3%. One might be tempted to attribute this dip 
to the occultation of a planet on an eccentric orbit, for which 
the occultation need not be halfway between transits. How¬ 
ever, the statistical significance of this dip is dubious, and the 
required value of the eccentricity would be e > 0.35, using 
Eq. 33 of Winn ( |2010| ). Such a high eccentricity would be 
unprecedented and~unexpected for a short-period planet. In 
general, giant planets with periods shorter than 3 days have 
nearly circular orbits, a fact that is attributed to the gradual 
action of tidal dissipation. Given the youth of the star, it is 
possible that there has not yet been sufficient time for orbital 
circularization; however, a higher eccentricity and a poten¬ 
tially smaller pericenter distance would also put the planet in 
even more danger of violating the Roche limit. 

Our non-detections of occultation signals at both 1.7 pm 
and 4.5 pm bands rule out the existence of a planet that radi¬ 
ates like a blackbody in these two bands. We have not pursued 
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t - t occ (hr) 

Fig. 7.— Top. —Near-infrared photometry of PTFO 8-8695 spanning the 
predicted time of occultation. The data points are averaged in groups of 5 
for clarity. The solid gray line shows the 3cr bound on the maximum depth 
of the occultation given by this light curve. The dashed grey line shows the 
minimum predicted occultation depth according to the model of Barnes et al. 
{2013}. Bottom. —Same, but after dividing through the light curve by the 
best-lit cubic function to the out-of-occultation region. 

more realistic models for the planetary emission spectrum, 
given that the atmospheric composition is unconstrained, but 
seems unlikely that atmospheric absorption features would 
suppress the planetary flux in both bands to such a degree that 
it would be undetectable in our data. 

2.4. Departure from periodicity 

Perhaps the most important finding of all the photometric 
observations is that the fading events are not strictly periodic. 
The top panel of Figure [TT] shows the residuals after subtract¬ 
ing the best-fitting linearTunction of epoch from the measured 
times of minimum light, 

t„ = t 0 +Pn , (4) 

for which x^ in = 15573 with 35 degrees of freedom. This poor 
fit is the result of the large scatter («15 min) of the residuals 
within each season, and the even larger deviation (^1.3 hours) 
of the most recent season’s residuals relative to the earlier 
data. The pattern of residuals suggests that the period was 
nearly constant up until the 2014/5 observing season, when 
the fading events began occurring earlier than expected. This 
apparent change in period or phase can be readily checked 
by gathering additional data over the next few seasons. The 
best-fitting parameters of the linear ephemeris are 

t 0 = 2455201.832 =b 0.007 days, 
p = 0.448391 ±0.000003 days. 

In these expressions the uncertainties have been scaled up by 
a factor of yj Xmin/^dof to account for the statistically poor fit. 

We also tried fitting a quadratic function of epoch, 

1 d P 9 

t n =t 0 +Pon+-—n 2 , (5) 

2 an 

for which x^ in = 4980 with 34 degrees of freedom. After 
enlarging the parameter uncertainties as described above, the 


best-fitting parameters are 

t 0 = 2455201.790 =b 0.006 days, 

P 0 = 0.448438 =b 0.000006 days, 
dP/dn = (-2.09 =b 0.25) x 10“ 8 days epoch -1 . 

The bottom panel of Figure [IT] shows the residuals between 
the observed and calculateamnes. The implied fractional 
change in period per epoch, calculated as is equal to 
-4.66 x 10 8 . If this period change were to continue steadily, 
the period would shrink to zero after Po/P ~ 10 4 years. 

3. SPECTROSCOPIC OBSERVATIONS 

We monitored the optical spectrum of PTFO 8-8695 on the 
night of 2013 Dec 12 UT, employing the Keck I telescope 
and its High Resolution Spectrograph (HIRES; Vogt et al. 
|1994| ). A total of 22 observations with integration times of 
14 min were obtained. This sequence covered 2 hr spanning 
the event, and 1 hr after the event. HIRES was used in its 
standard setting, but without an iodine cell in the light path. 

To confirm that a fading event was indeed taking place dur¬ 
ing the spectroscopic observations, we attempted to gather 
simultaneous photometry with several small telescopes, but 
in only one case was the weather at least somewhat coopera¬ 
tive. We obtained data in the H band with Mimir, a cryogenic, 
facility-class near-infrared instrument on the 1.83 m Perkins 
telescope outside Flagstaff, Arizona (Clemens et al. |2007| . 
Figure [12] shows the light curve. A transit-like dip of 
was seen at the expected time, confirming that a fading event 
did occur, although the data are too noisy to extract much fur¬ 
ther information. We also note that fading events were seen 
by TRAPPIST on 2013 Dec 11 and 14, bracketing our Keck 
observation. 

3.1. Search for the Rossiter-McLaughlin effect 

The primary purpose of the spectroscopic observations was 
to seek evidence for the Rossiter-McLaughlin (RM) effect, the 
spectroscopic anomaly that is seen during a planetary tran¬ 
sit due to stellar rotation. During a transit, a planet blocks 
different portions of the rotating stellar photosphere, leaving 
a particular imprint on the rotationally-broadened stellar ab¬ 
sorption lines. The exact shape and time development of the 
spectral deformations depend on the transit parameters, and in 
particular on the angle between the stellar rotation axis and the 
orbital axis as projected on the sky plane. For PTFO 8-8695, 
the planet hypothesis requires a large misalignment between 
these angles. We attempted to detect the RM effect and test 
that prediction. 

After the initial data reduction we corrected for the blaze 
function by using calibration lamp exposures to estimate the 
blaze function for each order, and then fitting a linear function 
of wavelength to remove the residual variations and normal¬ 
ize the continuum to unity. Following the barycentric correc¬ 
tion, all of the out-of-transit exposures were co-added to cre¬ 
ate a single spectrum with a higher signal-to-noise ratio. This 
was used for a final differential normalization, wherein the 
summed spectrum was subtracted from each observed spec¬ 
trum, and a 4th-order polynomial was fitted to the residuals in 
each order. These polynomials were subsequently subtracted 
from the corresponding spectrum. This was done to minimize 
the potential influence of any time variations in the blaze func¬ 
tion throughout the night. We verified that the details of this 
normalization process did not have a significant influence on 
the following analysis. 
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Fig. 8.— Left .—The top panel shows the expected occultation depth (5 OC c (solid black line) plotted as a function of the Bond albedo for / = 4. The red and 
blue lines represent the contributions to the total occultation depth from thermal emission of the planet and reflected light, respectively. For comparison, the 3a 
upper bounds on S oc c obtained from Magellan and Spitzer are shown as dashed and dotted lines. The bottom panel shows the corresponding dayside effective 
temperature of the planet. Right. —Same, but for / = 2. 
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Fig. 9.— Spitzer time-series 4.5 pm photometry of PTFO 8-8695 from 
2012 Apr 28. The small gray points represent individual measurements; the 
larger dark blue points are time averages. Red bars show the times of the 
fading event, and the predicted time of the planetary occultation. 



For each spectrum, we calculated the cross-correlation 
function (CCF) with reference to a synthetic spectrum. 
The synthetic spectrum was obtained from the PHOENIX 
database (Husser et al. |2013| ), for a star with T e ff = 3500 K, 
logg = 3.5 and solar metallicity. We selected the appropriate 
wavelength ranges for creating the CCF via a visual inspec¬ 
tion. We needed to locate areas for which the normalization 
seemed reliable, and where there were at least a few well- 
defined absorption lines. As had been reported by VE+12, 
there are only a few regions between 5000-7000 A suitable 
for this work. 

We calculated the mean CCF based on all of the out-of- 
transit observations. Then we subtracted the mean CCF from 
each individual CCF. When ordered in time, the resulting “dif¬ 
ferential CCFs” should display the shadow of the transiting 
object in velocity space. The deformation due to a transiting 
object would be seen as a dark line. The slope of this line in 
the velocity-time plane would depend on the projected obliq¬ 
uity. For example, in the case of good spin-orbit alignment, 
there would be a deficit of blue light (negative radial veloc¬ 
ities) in the first half of the transit, followed by a deficit of 
redshifted light during the second half. As can be seen in Fig¬ 


ure |T3j no such signal — neither aligned nor misaligned — 
can be discerned. 

Could we have detected the spectroscopic transit of the hy¬ 
pothetical planet, given the quality of our data? To answer this 
question, we simulated the RM signal of a transiting planet 
with (. R p /R *) 2 = 0.026, the approximate transit depth mea¬ 
sured by TRAPPIST during the events of 2013 Dec 11 and 
14, closely bracketing the event observed with Keck. We as¬ 
sumed vsin/* = 100 kms -1 , a value consistent with the line 
broadening seen in our Keck spectra, and adopted a macro- 


turbulent velocity of 15 kms -1 . Then we injected RM signals 
into the data, for various choices of A, the sky-projected stellar 
obliquity. 

Some representative examples of the simulated RM effect 
are shown in Figure [14] The left panel shows the simulated 
RM signal of a welFafigned planet (A = 0°). Such a signal 
would easily have been detectable with the data at hand. The 
same is true for A = 45°, shown in the middle panel. How¬ 
ever, for A = 90°, the signal would have been more difficult to 
detect. This is because in this case the RM signal is nearly sta¬ 
tionary in velocity, as shown in the right panel, making it more 
difficult to separate from the noise in the velocity-time plane. 
Based on visual inspection of figures similar to Figure [l4| we 
conclude that we can rule out any trajectory except fortnose 
within about 15° of perpendicularity (A = 90° or 270°). Given 
the non-Gaussian and correlated nature of the noise, it is dif¬ 
ficult to make a firmer statistical statement. 


3.2. Projected rotation rate 

A secondary goal was to seek changes in the sky-projected 
rotation rate (vsin/*) that would be expected if the star’s rota¬ 
tion axis is precessing around the total angular momentum of 
the system. If PTFO 8-8695 does consist of a star and a planet 
whose rotation axes and orbits precess around the common 
angular momentum, then v sin should change with time. We 
searched for such a change between the two epochs for which 
Keck/HIRES data have been obtained. Five spectra were ob¬ 
tained in April 2011 and presented by VE+12. Another epoch 
is represented by our December 2013 data. 

We derived CCFs from the 2011 data in exactly the same 
way as for the 2013 data. We compared the CCFs with theo¬ 
retical absorption lines taking into account uniform rotation, 
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Fig. 10.— Close-up of the “transit” and “occultation” data, after rectification to unit flux outside of the events. 
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Fig. 11.— Top panel : Timing residuals after subtracting the best-fitting 
linear function of epoch(constant period). Bottom panel : Timing residuals 
after subtracting the best-fitting quadratic function of epoch (steady decrease 
in period). 
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Fig. 13.— Time variations of the CCF during the night of 2013 Dec. 11/12. 
Dark areas indicate a deficit relative to the mean CCF on that night. The 
CCFs have been normalized to have a peak value of unity. The dashed lines 
indicate the predicted times of first and last contact. Vertical lines indicate 
our estimate for vsin /*. One can see variations of the CCFs before, during, 
and after the transit. No clear sign of a planetary transit is visible. 


limb darkening, macroturbulence, and gravity darkening. We 
adopted the quadratic limb darkening parameters from the ta¬ 
bles of Claret et al. p0l2| for logg = 3.5 and T eff = 3500 K, 
solar metallicity, and the Johnson V band. We adopted a 
macroturbulent velocity of 15 kms -1 . To model the gravity 
darkening we assumed the same effective temperature, a ro¬ 
tation period of 0.448 days, a stellar mass of 0.4 M 0 , a stel¬ 
lar radius of 1.4 R@, a stellar inclination angle of 90°, and 
a gravity-darkening exponent of /3 = 0.25. We neglected any 
oblateness of the stellar photosphere. As an example, Fig¬ 
ure |T5] shows the CCF for one of the pre-transit observations 
from Dec 2013, along with the best-fitting model. 

The results for vsin 4 are 103.6 ±0.3 kms -1 in 2011, and 
104.1 ±0.7 kms -1 in 2013. The quoted uncertainties are 
based on the scatter between the different observations for 
each epoch, and do not include any additional systematic un¬ 
certainties due to the limitations of the model (such as uncer¬ 
tainties in the treatment of limb darkening and gravity dark¬ 
ening, or the neglect of differential rotation and oblateness). 
Therefore the relative variation in vsin 4 is bounded to less 
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Fig. 14.— Same as Figure [13] but in this case an artificial signal of a transiting planet has been injected. The left panel shows the case A = 0° (spin-orbit 
alignment), the middle panel snows A = 45° and the right panel shows A = 90°. 



Fig. 15.— CCF for the pre-transit observation from Dec. 2013 (black), 
along with the best-fitting model with vsin/* = 103 kms -1 (blue). The red 
dashed line represents a model with vsim* = 80 kms -1 , the value reported 
by VE+12. The gray line shows the differences between the data and the 
best-fitting (blue) model, vertically offset by 0.3 units for clarity. 

than a percent, although the absolute value is probably uncer¬ 
tain by at least 10%. 

To obtain a better idea about the absolute value of v sin 4 
and its uncertainty, we tried fitting individual absorption lines 
rather than the CCF. Specifically, we fitted seven apparently 
isolated lines between 5300 and 7700 A. The standard de¬ 
viation in the vsin4 measurement form these seven absorp¬ 
tion lines varies between 7-15 kms -1 for the different ob¬ 
servations, which suggests that an uncertainty of 10 km s -1 in 
vsin4 should be a reasonable estimate. This leads to our final 
estimate of 103 =b 10 km s -1 . 

Our result for vsin4 is higher than the value of 80.6 db 
8.1 kms -1 reported by VE+12, but we do not think that this 
necessarily (or even likely) implies that the projected rota¬ 
tion rate is varying in time. This is because it is difficult to 
compare the results directly, given that VE+12 used a com¬ 
pletely different instrument and analysis technique. Our in¬ 
ternal comparison is much more sensitive, since it is between 
two Keck/HIRES spectra obtained at different times and ana¬ 
lyzed in exactly the same way. 

3.3. Time variations in emission-line profiles 


Our optical spectra also reveal strong emission lines from 
the hydrogen B aimer series as well as the Ca II H & K tran¬ 
sitions. The top panel of Figure [16] shows the median Ha 
line profile, based on all the “out-of-transit” spectra observed 
on 2013 Dec 12 (at least an hour before or after the time of 
minimum light). The line is very broad. Most of the emis¬ 
sion is confined to velocities <100 km s -1 , consistent with 
the star’s rotation rate, but the velocity profile extends to at 
least 300 km s -1 , particularly on the blue side. This is sug¬ 
gestive of at least a low level of ongoing accretion. Material 
that falls onto the star from large distances could attain the 


free-fall velocity yJlGM^/R* « 330 km s 1 , given the nomi¬ 
nal parameters M* = 0.4 M 0 and R* = 1.4 Rq. The equivalent 
width of the Ha line is 8.7 A, placing the star near the tradi¬ 
tional borderline between the categories of “weak-lined” and 
“classical” T Tauri stars. For simplicity it is often said that 
the classical stars are actively accreting while the weak-lined 
stars are not accreting, although in reality there seems to be no 
sharp distinction between these categories, and PTFO 8-8695 
presents an intermediate case. 

Figure [16] also shows the time sequence of observed 
changes inthe Ha line profile. Specifically, for each of the 
22 spectra, we plotted the residuals between the observed line 
profile and the median “out-of-transit” line profile. Also indi¬ 
cated are the times of minimum light, as well as “ingress” (one 
hour prior) and “egress” (one hour afterward). The red dashed 
line indicates a slope of 37 km s -1 hr -1 . This is the expected 
radial acceleration of any feature attached to the stellar photo¬ 
sphere, which would move from -100 km s -1 to +100 km s -1 
over the course of F rot /2 = 0.224 days. 

Evidently the line profile varied in a complex pattern on 
a timescale of minutes. Several excess-emission features do 
seem to be rotating along with the star; for example, a pattern 
of positive residuals appears in the fourth-to-final spectrum 
(time coordinate t = 1.015) at velocity +50 km s -1 and shifted 
redward at the expected rate throughout the final three ob¬ 
servations. A similar pattern - perhaps originating from the 
same feature - is seen starting at minimum light at velocity 
-50 km s _1 and proceeding redward until about an hour after 
minimum light. These particular components of the emission 
line seems likely to be caused by active regions on the stellar 
surface. 

During the fading event, the residuals show relative absorp¬ 
tion at a redshifted velocity of 25-100 km s -1 . The absorption 
seemed to disappear at around the same time as the end of 
the fading event. The onset of the absorption feature was at 
least 2 hours before minimum light, which is at least an hour 
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before what seems to be “ingress” of the fading event. Thus, 
the transient redshifted absorption does not seem to be exactly 
coincident with the fading event, although it does at least sug¬ 
gest that the fading event was associated with hot infalling 
material in front of the star. 

The H 7 line profiles (not shown here) tell a similar story 
but with a lower signal-to-noise ratio. The H [3 line was not 
observed, given the spectral format. Figure [TT] shows the me¬ 
dian Ca II H line profile, along with the time series of devi¬ 
ations from the median. In this case the fractional variations 
were even stronger and seemingly faster; there is no straight¬ 
forward narrative to the sequence of residuals. 

To seek independent evidence for ongoing accretion, we 
used the available broadband photometry to construct the 
spectral energy distribution of PTFO 8-8695 and search for 
any infrared excess. Figure |T8] shows the result, based on a 
query of the VizieR websitq^j which gave measurements in 
the VRJHK bands as well as the WISE W1-W4 bands. (The 
W 4 observation gave an upper limit.) We corrected for dust 
extinction with the dust map from the NAS A/IPAC websit^J 
and fitted the results to a grid of zero-metallicity stellar atmo¬ 
sphere models from the library of Castelli & Kurucz (|2004|. 
The best-fitting stellar parameters were r eff = 3500 ±120 K 
and logg = 4.0 ± 0.9. The effective temperature is in agree¬ 
ment with the previously reported value of 3470 K (Briceno 
et al. |2005j ). The apparent lack of an infrared excess out to 
10 fim islmaracteristic of weak-lined T Tauri stars. This lack 
of evidence for the existence of an accretion accretion disk 
within 1 AU does not necessarily rule out accretion, but does 
suggest that any accretion is relatively weak. 


4. DISCUSSION 

We now summarize the results of the three tests that we 
undertook to test the hypothesis that the fading events of 
PTFO 8-8695 are the transits of a close-in giant planet: 

1. The new light curves show variations in depth and 
duration from event to event. However we did not 
find strong evidence for asymmetries or other secular 
changes in morphology indicative of the changing tra¬ 
jectory of a transiting planet. Furthermore, in all cases, 
a fading event was observed at the appointed time, even 
though the cessation of transits was predicted to be 
likely by Barnes et al. ( |2013| ). 

2. Infrared photometry spanning the predicted times of 
planetary occultations has ruled out signals of the ex¬ 
pected amplitude. 

3. The Rossiter-McLaughlin effect could not be detected, 
ruling out a transiting planet with the expected param¬ 
eters, unless the planet’s trajectory is nearly perpendic¬ 
ular to the projected stellar equator. Nor did we detect 
any changes in vsin± between 2011 and 2013, which 
would have been produced by precession of the stellar 
rotation axis. 


Any of these tests could have resulted in a confirmation of 
the planet hypothesis. In all cases, though, the planet hypothe¬ 
sis was disfavored. In addition we found that the fading events 
are not strictly periodic. In the most recent observing season 
the events occurred more than one hour earlier than expected. 

18 http://vizier.u-strasbg.fr 

19 http://irsa.ipac.caltech.edu 
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Fig. 16.— The Ha line profile of PTFO 8-8695 on 2013 Dec 12. Top.— 
Median line profile, based on all spectra more than one hour away from min¬ 
imum light. Bottom .—Time series of residuals between each observed spec¬ 
trum and the median. Open symbols are the “out-of-transit” spectra used to 
create the median spectrum. Filled symbols are the spectra within one hour 
of minimum light. On the left axis, each vertical tick mark represents one 
unit of relative flux, on the same scale as the top panel. The right axis gives 
the time of each spectrum. 


This finding is incompatible with the strict periodicity that one 
expects for a planetary orbit. 

While this paper was in preparation we learned of the work 
by Ciardi et al. ( |2015| ) and Koen ( [2015) ), who also pursued 
follow-up observations of PTFO 8-8695. Among the observa¬ 
tions of Ciardi et al. ( |2015| ) was a nondetection of any transit¬ 
like event on 2012 Dec 21, based on observations in the r' 
band. They gave an upper limit of 0.7% on the transit depth. 
Their simultaneous spectroscopy also revealed no evidence of 
the Rossiter-McLaughlin effect. They interpreted these non- 
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Fig. 17.— TheCaIIHlineprofileofPTFO8-8695 on2013Dec 12. Top.— 
Median line profile, based on all spectra more than one hour away from min¬ 
imum light. Bottom .—Time series of residuals between each observed spec¬ 
trum and the median. Open symbols are the “out-of-transit” spectra used to 
create the median spectrum. Filled symbols are the spectra within one hour 
of minimum light. On the left axis, each vertical tick mark represents 5 units 
of relative flux, on the same scale as the top panel. The right axis gives the 
time of each spectrum. 


detections as evidence for the predicted cessation of the tran¬ 
sits. However, our observations reveal that the fading events 
did indeed take place on 2012 Dec 11,14, and 15, with a depth 
of approximately 3% in all cases. It seems unlikely that or¬ 
bital precession could have abruptly reduced the transit depth 
from 3% to below 0.7% in less than one week. Therefore our 
results cast doubt on this aspect of th e interpretation of Ciardi 
et al. ( |2015b . Likewise, Koen ( 2015| ) reported non-detections 
of predicted fading events on 2015 Jan 3-4. However, the pre¬ 
dicted times were based on the assumption of strict periodic- 



Fig. 18.— Spectral energy distribution of PTFO 8-8695 based on pub¬ 
licly available broadband photometry (orange diamonds), along with the best¬ 
fitting stellar atmosphere model (blue line). The data were corrected for ex¬ 
tinction (red diamonds) prior to fitting. There is no evidence for any infrared 
excess that could be attributed to a circumstellar disk. 


ity, which our observations have shown to be false. Judging 
from Fig. 2 of Koen ( |2015| ), it seems possible that the fading 
events were recorded in the SAAO observations a few hours 
earlier than expected. In any case we detected a clear 2% 
fading event on 2014 Dec 27, only one week earlier than the 
SAAO observations. 

In summary, our observations have significantly reduced 
the credibility of the planet hypothesis. However, because 
the hypothesis invokes an unusual planet in unusual circum¬ 
stances, it is difficult to rule out definitively. It may be pos¬ 
sible to find reasons for the failure of each of the individual 
tests. For example, the predicted occultation times might have 
been inc orrect, because the planet has a highly eccentric or¬ 
bit (see § |2.3.2| ). Or perhaps the values of the key parameters 
( R p /R *) 2 and R*/a are smaller than the values postulated by 
Barnes et al. (1201 3j ), which would reduce the predicted oc¬ 
cultation depth. The planet’s atmosphere might have deep 
absorption features near 1.7 fim and 4.5 fi m that rendered 
the planetary occultations undetectable. The planet’s orbit 
might have been nearly perpendicular to the stellar equator 
at the time of our attempt to detect the Rossiter-McLaughlin 
effect. It remains possible that a comprehensive search of pa¬ 
rameter space of the model proposed by Barnes et al. ( 2013 ) 
— including the effects of gravity darkening and orbital pre¬ 
cession — could reveal a configuration that possesses these 
properties and is also compatible with the lack of detectable 
change in v sin 4 , as well as the unexpectedly bland morpholo¬ 
gies, strongly chromatic depths, lack of occultation signals, 
and timing irregularities that are seen in the new light curves. 
We leave such a computationally intensive search for future 
work. It is also important to try and develop alternative hy¬ 
potheses for the fading events of PTFO 8-8695. Below we 
describe four alternatives, along with their apparent strengths 
and weaknesses. Figure [19] illustrates these hypotheses. 


4.1. Dust-emitting rocky planet 

The seemingly erratic variations in the depth and duration 
of the fading events, along with the slight asymmetries in the 
phase-folded light curves (Figure^, bring to mind the case of 
KIC 12557548 (Rappaport et al. ||2012( ). This object was iden¬ 
tified in data from the Kepler spacecraft. It exhibits transit¬ 
like flux dips with a very short period (15.7 hours) and du¬ 
ration (1.5 hours), an erratically varying depth ranging from 
0-1.3%, and an egress of longer duration than ingress. Rap- 
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Fig. 19.— Illustrations of the five hypotheses discussed in Section]?] 

paport et al. ( [2012| interpreted the dips as transits by a dusty 
tail being emitted by a small rocky planet. Two other sim¬ 
ilar cases have since been identified (Rappaport et al. [2014} 
Sanchis-Ojeda et al. |2015| ). In at least one case the transit 
depth has been shown to be smaller at infrared wavelengths, 
as expected for extinction by small dust grains. 

PTFO 8-8695 shares much of the phenomenology that has 
just been described. The main difference is that the other sys¬ 
tems are not rapidly-rotating young stars. They appear to be 
slowly-rotating main-sequence stars, and are not even close 
to being synchronized with the transit period. Furthermore it 
is not clear whether this hypothesis could be reconciled with 
the apparent change in period that was seen in the most recent 
season of observations. 

It is tempting to invoke tidal dissipation as a mechanism for 
gradually shrinking the orbit and shortening the period, but 
this is implausible for a low-mass rocky body. Conceivably, 
orbital decay could be a consequence of the dust emission. If 
the dust acquires additional specific angular momentum from 
radiation pressure while leaving the system, it would be driven 
into a higher orbit. There it would pull back on the planet 
and potentially extract angular momentum from the planetary 
orbit. However, the magnitude of this effective drag is diffi¬ 
cult to estimate from first principles, particularly because the 
dust may represent only a modest fraction of the total mass 
loss, and the gas need not behave the same way as dust as it 
leaves the system. In any case t he lif etime of ~10 4 years im¬ 
plied by the observations (see § |2.4| ) is uncomfortably short. 
It would require a special coincidence to observe such a short¬ 
lived phase of evolution. 


4.2. Starspots 

The synchronization of the stellar rotation period and the 
period of the fading events, along with the changing depth 
and duration of the fading events, raises the suspicion that the 
periodic dips are caused by starspots being carried around by 
rotation. The star is expected to be heavily spotted, given its 
youth. Moreover, the photometric variations produced by stel¬ 


Accretion hotspot 



lar activity are expected to be weaker in the infrared than in 
the optical, consistent with our observations. Gradual changes 
in the spot pattern could be invoked to explain the changes in 
depth, duration, and timing of the fading events. 

VE+12 have already pointed out the main weakness of this 
hypothesis. Flux variations caused by starspots have a nat¬ 
ural timescale of half the rotation period, the interval over 
which a patch of the stellar surface is typically visible to the 
observer. In contrast, the fading events of PTFO 8-8695 last 
only «15% of the rotation period. One can arrange for short- 
duration dips by locating the starspots near one of the rotation 
poles, and tipping the star such that the spots are only on the 
visible hemisphere over a narrow range of longitudes. Indeed, 
T Tauri stars are well known for displaying stable spot pat¬ 
terns near the magnetic poles (e.g. Granzer et al. 2000| . How¬ 
ever, in such a circumstance the spot is highly foreshortened 
and limb-darkened, and it is difficult to produce variations of 
several percent. It is also difficult to produce the apparently 
sharp features that have been occasionally seen in the light 
curves, such as the apparent point of “fourth contact” in our 
Magellan if-band time series (see the bottom panel of Fig¬ 
ure [6]). It might be possible to save this hypothesis by invok¬ 
ing multiple spots and in a complex and stable pattern, such 
that the summation of their photometric variations is coinci¬ 
dentally narrow in time. This model, though, would be rather 
contrived. 


4.3. Eclipses by a circumstellar disk or dust 

As we have just seen, the basic difficulty of any model in 
which the variations are due to features on the stellar pho¬ 
tosphere is the relatively short duration of the fading events. 
Eclipses by orbiting bodies avoid this problem because the 
time spent in front of the star scales as R*/a, and the orbital 
distance a can be adjusted to match the observations. This 
is the basis of the planet hypothesis, for which Barnes et al. 
(2013) found R*/a = 0.5-0.6. However the eclipsing body 
need not be a planet. Perhaps it is a feature in the innermost, 
corotating portion of the stellar magnetosphere. 
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At an age of a few million years, low-mass T Tauri stars are 
often still actively accreting material from their circumstellar 
disks. PTFO 8-8 695 l acks any evidence for an accretion disk, 
yet, as noted in § 33 the strength and breadth of the Ha line 
profile places the star in between the traditional categories of 
“classical” and “weak-lined” T Tauri stars and suggest that 
PTFO 8-8695 may be weakly accreting. So we will consider 
the implications of a small amount of dust that would not be 
detected as an IR excess, but would still be able to produce 
the flux dips. The accretion process is thought to proceed as 
follows (see, e.g., Bouvier et al. [2007| for a review). First, 
matter spirals inward through a thin disk, until it reaches an 
orbital distance of a few stellar radii, at which point the disk is 
disrupted by the star’s magnetosphere. Then, within the mag¬ 
netosphere, the material is magnetically funneled into narrow 
tubes or columns, extending from the inner edge of the ac¬ 
cretion disk onto the star’s magnetic poles. The material falls 
freely along these columns and crashes onto the star, produc¬ 
ing shock waves and a luminous hotspot. In many models, the 
stellar rotation rate becomes synchronized with the Keplerian 
orbital velocity at the inner radius of the magnetosphere (the 
corotation radius), a phenomenon known as disk locking. 

This picture contains several elements that could naturally 
lead to flux variations with a period equal to the stellar rota¬ 
tion period. For example if the accretion disk is warped or has 
other non-axisymmetric variations in thickness near the inner¬ 
most, corotating portion of the disk, then the star may be peri¬ 
odically eclipsed by these irregularities. This is thought to be 
the basic explanation for the quasiperiodic eclipses of AA Tau 
(Bouvier et al. |1999| 20031 ). Alternatively, periodic eclipses 
could be produced by stable patterns or concentrations of dust 
in the accretion flow. This may explain the observed variabil¬ 
ity of the “short-duration dippers” recently identified by Cody 
et al. ( |2014| ) and Stauffer et al. ( |2015| . Those authors studied 
a number of T Tauri stars in the young cluster NGC 2264 
that exhibit short-duration fading events. They found that the 
flux dips are quasi-periodic and exhibit changes in depth and 
shape from epoch to epoch over a period of several years. The 
dips appear shallower in the infrared than in the optical, and 
the light curves have rounded minima rather than being flat- 
bottomed (Stauffer et al. 2015| ). The flux dip periods are usu¬ 
ally equal to the rotation periods of the stars, and some dips 
have been observed to persist for years. All these properties 
have been attributed to extinction by infalling material from 
the innermost portion of the accretion disk (McGinnis et al. 
2015| ). And, all these properties are consistent with our obser¬ 
vations of PTFO 8-8695. 

In some respects, though, PTFO 8-8695 is different from 
the rest of the dippers. Its rotation period of 0.45 days is 
shorter than the 3-10 day pe riods of most of the stars ob¬ 
served by Cody et al. ( |2014| ) and Stauffer et al. ( 2015[ ), or 
the 8 day period of A A Tau. The duration of the fading events 
is also relatively shorter at 15% of the period, compared to 
the more typical value of ~30%. It is also unclear whether 
disk warps or dust concentrations could produce the appar¬ 
ently sharp features and flat bottoms observed in a few of our 
light curves. Furthermore, the dippers all have SEDs with a 
detectable infrared excess, but the currently available data for 
PTFO 8-8695 show no evidence for any infrared excess (Fig¬ 
ure [18]). 

Furthermore, it is questionable whether dust can exist in 
solid form so close to the star, with an orbital distance less 
than 2 R+. Assuming that stellar radiation is the dominant 
mechanism of heating, the dust sublimation radius R s is given 


by (Monnier & Millan-Gabet |2002| 


Rs 2 VQr f T 


R * 


( 6 ) 


where T s is the dust sublimation temperature (^1500 K), and 
Qr = Gabs(?*)/Gabs (T s ) is the ratio of the dust absorption effi¬ 
ciencies for radiation at the color temperatures of the incident 
and reemitted fields, respectively. Adopting the customary 
value of T s = 1500 K, and assuming Q R to be of order unity 
(as expected in this case for silicate grains), this formula gives 
R s = 2.7 R *. Therefore the hypothesized dust with R*/a= 0.5- 
0.6, corresponding to a =1.7-2 R+, would likely be vaporized. 

4.4. Occultations of an accretion hotspot 

Hotspots are another aspect of the magnetospheric accre¬ 
tion model that has previously been invoked to explain some 
of the quasiperiodic variations of T Tauri stars (Herbst et al. 
|1994| ). We advance here a related hypothesis for the case of 
PTFO 8-8695: perhaps the fading events represent occulta¬ 
tions of one of the hotspots that is produced by ongoing low- 
level accretion. In this scenario, the accreting material is fun¬ 
nelled onto a hotspot near one of the star’s magnetic poles, 
which is displaced from the star’s rotation pole. Furthermore, 
the star’s rotation pole is tipped toward the observer such that 
the hotspot is on the visible hemisphere for «85% of each 
rotation period. When the hotspot is hidden from view, we 
observe a fading event. 


This is similar to the starspot hypothesis (§ 4.2) in that the 
photometric modulations are the result of the rotation of a 
photospheric feature of contrasting intensity. However, re¬ 
placing the dark starspot with a luminous hotspot could po¬ 
tentially solve some of the problems faced by the starspot 
model. Hotspots have been inferred to have effective temper¬ 
atures ranging up to 10 4 K, and thereby present much higher 
contrast than starspots with the surrounding photosphere. Fur¬ 
thermore, hotspots are probably not confined to a vertically 
thin surface layer of the photosphere, and as such they are 
not subject to the effects of limb-darkening or foreshortening. 
Together these factors may make it easier to produce abrupt 
modulations of a few percent in the total light even for a small 
feature near the stellar limb. 

If this model is correct, then the mass accretion rate can be 
estimated from the loss of light during fading events, which 
is observed to be of order 5%. This requires the accretion 
luminosity to be 


GMM 

~~R 


' 0.05 L bol . 


Using the nominal values M * = 0.4 M 0 , R = 1.4 Rq, and 
L bo i = 0.25 L e , we find M ~ 10 -9 Mq yr -1 , at the low end 
of the range of mass accretion rates that has been inferred for 
accreting T Tauri stars (10 -9 -10 -7 M 0 yr -1 ). This seems rea¬ 
sonable: a relatively low accretion rate is also consistent with 
the relatively weak Ha equivalent width of 8.7 A and the ab¬ 
sence of a detectable infrared excess. 

One concern with this model is that in the systems where 
hotspots have been previously invoked to explain photomet¬ 
ric variability, the pattern of variability is not as consistent or 
long-lasting as is seen in PTFO 8-8695 (Herbst et al. |1994| . 
Oftentimes the hotspot variability shows no periodicityor, at 
best, short-lived periodicity, sometimes with period changes 
of 20% or more. However, those previous results pertain to 
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classical T Tauri stars with much higher inferred accretion 
rates; perhaps we are seeing the different and more stable be¬ 
havior of a hotspot in a more weakly accreting system. 

4.5. Summary 

We have discussed five hypotheses for the fading events of 
PTFO 8-8695. The precessing giant-planet model has failed 
several key tests, the most serious of which are probably the 
nondetection of the planetary occultations, and the apparent 
shift in the phase of the fading events in the most recent sea¬ 
son of observations. The planet hypothesis also struggles to 
explain the observed coincidence of the rotation and orbital 
periods. This same problem afflicts the hypothesis of the dust- 
emitting rocky planet. 

The other three models share the virtue of a natural expla¬ 
nation for the equality of the rotation period and the period 
of the fading events. However, the starspot model has dif¬ 
ficulty reproducing the observed duration and occasionally 
sharp ingress/egress of the fading events. The other two mod¬ 
els invoke the presence of an accretion disk, for which the 
evidence is ambiguous or negative: the Ha line profile does 
extend to higher velocities than can be explained by stellar ro¬ 
tation, but the equivalent width is relatively low and there is 
no detectable infrared excess. 

The occulted-hotspot model seems quantitatively promis¬ 
ing, as it is consistent with a low accretion rate of 
^10 -9 Mq yr -1 . There is no deterministic theory for the ex¬ 
pected photometric variations due to magnetically-funnelled 
accretion, making it difficult to achieve a firm confirmation 
of this hypothesis. Nevertheless all our observations seem at 
least consistent with this picture. At present, we consider this 
hypothesis to be the best explanation for PTFO 8-8695. To 
come to a firmer conclusion will probably require more pho¬ 
tometric and spectroscopic observations, seeking changes in 
the timing and behavior of the fading events, variations in the 
Ha line profile, and more sensitive searches for any infrared 
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(or ultraviolet) excess or other indicators of low-level accre¬ 
tion. 

At the outset of this project, and throughout most of this 
paper, we have been mainly concerned with the status of the 
planetary hypothesis for this intriguing planetary candidate. 
In fact this object may turn out to be useful for understand¬ 
ing magnetospheric accretion, due to a fortuitous geometry, 
thereby joining the ranks of such systems as A A Tau (Bou¬ 
vier et al. |1999| and KH 15D (Hamilton et al. [2012| . 
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